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ABSTRACT 


Continental ity:  Its  Estimation  and  Physical 
Significance.  (December  1985) 

Juan  Manuel  Yee  Fong,  B.S.,  University  of  New  Mexico 
.  Chairman  of  Advisory  Committee:  Dr.  Dennis  M.  Driscoll 

An  attempt  has  been  made  to  quantify  continental ity  better  than 
exists  in  current  indexes.  Conventional  indexes  of  continental ity  were 
examined  and  their  deficiencies  noted.  Several  different  approaches 
were  taken  to  develop  alternative  continental ity  indexes.  These  ap¬ 
proaches  were:  dividing  annual  temperature  range  by  the  difference  be¬ 
tween  summer  and  winter  insolation,  regressing  annual  temperature  range 
on  latitude  and  isolating  the  residuals  such  that  positive  ones  repre¬ 
sent  continental  locations  $pd  negative  ones  indicate  maritime,  or  less 
continental,  locations;  and  finding  the  summer  and  winter  lag  of  temper¬ 
ature  behind  radiation  with  the  cubic  spline  interpolation  technique. 

All  three  approaches  were  used  to  examine  the  North  American  continent 
as  a  whole.  In  addition,  the  Rocky  Mountain  and  Great  Lakes-Appalachia 
regions  were  chosen  for  mesoscale  analysis  in  order  to  determine  what 
effect  elevation  and  large  water  bodies  have  on  lag. 

Isopleth  analyses  of  the  plotted  values  of  the  proposed  alternative 
indexes  were  compared  with  conventional  index  patterns  for  North  America 
and  with  each  other.  Each  newly  developed  index  was  evaluated  as  an 
appropriate  measure  of  continental ity  and  a  physical  (meteorological) 
explanation  attempted  for  the  patterns. 
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1.  INTRODUCTION 

Continental ity  is  a  comparative  measure  of  the  differences  between 
climates  controlled  by  land  masses  and  those  controlled  by  water.  It 
was  first  used,  apparently,  by  W.  Zenker  in  the  1890' s.  Continental ity 
reflects  the  influence  of  the  strength  and  direction  of  the  general 
atmospheric  circulation,  topography,  and  aridity  (Griffiths  and 
Driscoll,  1982).  The  degree  of  continental ity  of  a  given  location  is, 
in  general,  proportional  to  its  distance  from  a  large  body  of  water, 
modified  by  prevailing  wind  direction  and  topography.  Oceanic ity  is  the 
opposite  of  continental ity,  and  characterizes  a  marine  climate. 

As  a  concept,  continental ity  can  be  useful  in  helping  to  explain 
the  similarities  and  differences  in  climate  between  various  locations. 
There  is  an  inherent  problem,  however,  in  defining  continental ity  in  a 
detailed  and  unambiguous  manner.  Part  of  the  problem  is  that  there  is 
no  way  of  assessing  absolute  continental ity. 

There  are  several  ways  in  which  continental ity  has  been  defined. 

One  method  is  from  a  geographic  point  of  view.  The  closer  a  location  is 
to  the  center  of  a  land  mass,  the  more  continental  it  is.  Another  defi¬ 
nition  can  be  based  on  a  location's  distance  from  the  nearest  ocean  or 
other  large  body  of  water.  A  third  way  of  describing  continental ity  is 
in  terms  of  climatic  extremes  and  variability  (e.g.,  temperature).  A 
continental  location  experiences  more  climatic  extremes  and  variability, 
and  has  more  seasons  than  oceanic  locations,  where  there  is  often  only 
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one  season  due  to  the  uniformity  of  the  climate.  These  definitions  are 
by  their  nature  relative  measures  of  continental ity,  since  they  stress 
one  or  more  variables,  no  one  of  which  gives  the  absolute  continental ity 
index. 

Of  all  climatic  elements,  temperature  is  probably  the  one  most  pro¬ 
foundly  influenced  by  continental ity.  The  most  conspicuous  difference 
between  regions  having  continental  climates  and  those  with  maritime 
climates  is  reflected  in  their  daily,  seasonal,  and  annual  temperature 
ranges  as  well  as  in  their  lag  times  between  radiation  and  temperature. 
Other  meteorological  elements  affected  by  continental ity  are  vapor 
pressure,  cloud  cover,  precipitation,  and  wind  direction  and  speed. 

There  is  more  evaporation  over  water  than  over  land  and  therefore  humi¬ 
dities  are  lower  in  continental  climates.  With  less  evaporation  comes 
reduced  cloudiness  and  rainfall.  Air  masses  originating  over  continen¬ 
tal  interiors  form  high  pressure  systems  around  which  wind  flow  is  anti- 
cyclonic  (clockwise).  Greater  surface  roughness  makes  wind  speeds  lower 
over  land  areas  thar  over  relatively  smooth  oceans. 

The  physical  basis  of  continental  ity  "is  in  the  different  thermal 
properties  of  water  and  the  various  land  substances.  Soils  have  low 
heat  capacities  and  heat  transfer  is  through  molecular  conduction.  In¬ 
coming  solar  radiation  (insolation)  penetrates  only  a  very  shallow 
layer.  Virtually  all  heat  energy,  therefore,  is  gained  and  lost  at  the 
surface.  For  this  reason,  land  areas  heat  and  cool  quickly.  The  air 
above,  having  a  very  low  heat  capacity  and  very  large  internal  mobility, 
also  changes  its  temperature  rapidly  in  response  to  the  heating  and 
cooling  at  the  surface.  Water,  on  the  other  hand,  experiences  smaller 


temperature  fluctuations  due  to  its  high  heat  capacity  and  its  ability 
to  transport  heat  via  convection.  Furthermore,  insolation  penetrates 
through  a  much  deeper  layer  than  in  soil;  not  very  much  energy  remains 
to  heat  the  air  above.  In  addition,  evaporation  is  greater  over  water 
than  over  land  surfaces,  leaving  less  insolation  available  to  heat  the 
water.  Lower  humidities  contribute  to  larger  temperature  ranges  over 
land.  As  a  consequence,  daily  and  annual  temperature  ranges  are  greater 
over  land  than  over  water. 
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2.  BACKGROUND  AND  LITERATURE  REVIEW 

Quantitative,  or  semi -quantitative,  measures  of  continental ity  have 
been  based  on  characteristics  of  radiation,  temperature,  precipitation, 
and  frequency  of  air  mass  types.  Continental ity  is  manifested  in 
temperature  in  three  principal  ways:  range,  lag,  and  variability. 

Range  and  lag  have  been  used  in  the  construction  of  continental ity  in¬ 
dexes.  Indexes  based  on  range  have  the  general  form 

k  =  (mA/(sin0  +  b))  +  n,  (1) 

where  k  is  the  index  of  continental ity,  A  the  average  annual  temperature 
range  (difference  between  the  mean  monthly  temperatures  of  the  extreme 
months,  usually  July  and  January  in  extra-tropical  areas),  0  the  lati¬ 
tude;  and  b,  m,  and  n  are  constants.  Indexes  based  on  air  masses  re¬ 
flect,  of  course,  the  frequency  of  occurrence  of  air  masses  originating 
over  continents,  as  opposed  to  those  originating  over  large  bodies  of 
water. 

a.  Previous  studies 

Brunt  (1924),  following  studies  of  R.  Spitaler  and  G.  Swoboda,  re¬ 
lated  continental ity  to  the  annual  temperature  range  and  the  annual 
range  of  the  average  global  solar  radiation  intensity.  His  equation  is: 

Ne  +  0.12  =  At/ (1 30.61  AS) ,  (2) 

where  Ne  is  the  continental ity  index.  At  is  the  annual  range  of  mean 
monthly  temperature  (July  minus  January),  and  AS  is  the  difference  be¬ 
tween  the  mean  global  solar  radiation  accumulated  in  the  months  of 


January  and  July.  Curiously,  Brunt  did  not  include  a  map  of  his  re¬ 
sults,  nor  did  he  specify  units  for  temperature  or  radiation.  The  quan¬ 
tity  Ne  +  0.12  is  equivalent  to  the  ratio  of  the  observed  annual  temper¬ 
ature  range  to  a  standard  range  of  temperature  expected  for  the  loca¬ 
tion’s  latitude.  It  expresses  the  fraction  of  a  1°  latitude  band  that 
would  have  to  be  covered  by  land  in  order  for  the  entire  parallel  of 
latitude  to  have  the  same  annual  temperature  range  as  that  observed  at  a 
given  latitude.  For  example,  a  station  at  45°N  and  Ne  =  0.5  means  that 
half  of  the  45°  parallel  must  be  covered  by  land  in  order  that  the  mean 
annual  temperature  range  of  the  entire  parallel  be  the  same  as  the 
station's  observed  annual  range.  This  term  also  measures  the  response 
in  monthly  mean  temperature  to  a  unit  change  in  the  mean  solar  radiation 
intensity.  If  the  response  is  zero  (e.g.,  an  oceanic  climate),  then  the 
corresponding  continental ity  index  is  -0.12.  Calculated  values  of  Ne 
ranged  from  about  -0.12  in  mid-ocean  in  the  southern  hemisphere  to 
approximately  1.3  in  Siberia. 

Zenker  (Conrad,  1946)  proposed  the  following  formula  for  conti¬ 
nental  ity: 


k  =  (1.2A/0)  -  20  .  (3) 

Latitude  was  included  to  account  for  generally  increasing  annual  temper¬ 
ature  range  with  increasing  latitude.  Johansson,  Gorczynski,  and  others 
modified  Zenker's  formula  by  taking  the  sine  of  latitude  and  introducing 
different  constants.  The  constants  for  Johansson's  formula. 


k  =  (1 .6A/sin0)  -  14 


(4) 


were  derived  by  arbitrarily  choosing  two  stations  which  seem  to  repre¬ 
sent  best  a  totally  continental  and  a  totally  oceanic  climate 
(Landsberg,  1958).  Stations  selected  were  Thornshavn  in  the  Faroe 
Islands  (k  =  0%)  and  Verkhoyansk  in  eastern  Siberia  (k  =  100%).  At  low 
latitudes,  k  may  far  exceed  its  theoretical  maximum  of  100%.  In  some 
cases  k  is  negative.  Temperature  ranges  are  so  small  in  the  interval 
1 0°N  -  10°S  that  latitude  in  this  equation  becomes  the  dominant  factor 
in  the  expression  for  continental ity. 

Conrad  (1946)  introduced  a  correction  factor  that  yields  more 
reasonable  values  for  the  tropics.  Ten  degrees  are  added  to  the  lati¬ 
tude  to  give  the  following  equation: 

k  =  (1.7A/sin(0  +  10))  -  14  .  (5) 

Conrad’s  formula,  however,  still  produces  some  negative  values.  Another 
disadvantage  is  that  it  loses  its  validity  at  latitudes  greater  than  80 
degrees.  To  correct  for  this,  for  all  latitudes  greater  than  80  de¬ 
grees,  (0  +  10)  is  assumed  to  equal  90. 

Several  workers  have  plotted  and  analyzed  continental ity  according 
to  Conrad's  definition  for  a  number  of  areas.  These  areas  include:  New 
England  (Fobes,  1954),  the  western  United  States  (D'Ooge,  1955),  the 
Great  Lakes  region  (Kopec,  1965),  the  Texas  coastal  zone  (Jehn,  1977a 
and  1977b),  and  Anglo-America  (Trewartha,  1981).  Using  Johansson's  for¬ 
mula  (Eq.  4),  MacKay  and  Cook  (1963)  constructed  a  map  of  continental ity 
for  Canada.  Each  of  these  maps  shows  the  moderating  influence  of  upwind 
bodies  of  water  and  distance  from  that  water.  Effects  of  topography 
also  are  apparent.  In  North  America,  index  values  are  generally  lower 


east  of  the  Great  Lakes  and  Hudson  Bay  and  are  higher  to  the  west.  Iso¬ 
lines  of  continental ity  are  more  closely  spaced  along  the  coastlines  of 
North  America  than  in  the  interior.  In  addition,  the  gradient  along  the 
west  coast  of  the  continent  is  tighter  than  along  the  eastern  seaboard 
(Fig.  1). 

Prescott  and  Collins  (1951)  constructed  worldwide  maps  of  the  lag 
of  air  and  sea  surface  temperatures  behind  radiation  by  means  of  elemen¬ 
tary  Fourier  analysis.  They  calculated  the  phase  angle  of  the  first 
harmonic  of  the  annual  temperature  curve,  and  that  of  twelve  months  of 
solar  radiation,  and  expressed  the  difference  in  degrees,  which  is 
essentially  the  same  as  days.  This  lag  of  temperature  behind  solar 
radiation  in  days  is  in  itself  an  index  of  continental ity.  Coastal  sta¬ 
tions  are  late  in  phase  while  continental  stations  are  early.  Lags 
along  the  western  edge  of  continents  in  the  middle  latitudes  are  longer 
than  along  the  eastern  edges,  due  to  the  influence  of  the  westerly  flow. 

In  a  continental  climate,  spring  months  and  the  corresponding  autum 
months  have  about  the  same  mean  temperatures.  At  coastal  stations, 
months  with  comparable  temperatures  occur  later  in  the  year.  This 
delay,  or  lag,  in  maritime  climates  was  used  by  F.  Kerner  to  define  an 
index  of  oceanicity  called  the  thermoisodromic  ratio: 

0  =  ((To  -  Ta)/A)100  (6) 

where  0  is  the  degree  of  oceanicity  in  percent,  A  is  the  annual  temper¬ 
ature  range,  and  To  and  Ta  are  the  October  and  April  mean  temperatures, 
respectively  (Landsberg,  1958).  Kerner's  results  agree  fairly  well  with 
conventional  maps  of  continental ity.  A  high  degree  of  oceanicity 
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corresponds  to  low  continental ity  index  values. 

Sumner  (1953)  calculated  the  standard  deviation  of  mean  monthly 
temperatures  for  a  number  of  stations  in  Anglo-America.  He  did  not  use 
temperature  variation  as  an  index  of  continental ity,  but  his  patterns  of 
mean  monthly  temperature  standard  deviations  are  very  similar  to  those 
shown  by,  for  example,  Conrad's  map  of  continental ity  for  North  America. 

Oliver  (1970)  plotted  average  temperature  and  precipitation  data 
for  all  twelve  months  for  a  number  of  stations  on  a  climograph  and  con¬ 
nected  the  extreme  points  with  a  straight  line  to  obtain  the  long  axis. 
The  longer  the  axis,  the  larger  is  the  annual  temperature  range  and  the 
more  continental  is  the  station.  Also,  a  small  angle  (subtended  from 
the  vertical)  indicates  the  dominance  of  continental  air  masses  (those 
originating  in  the  interior  of  a  continent)  throughout  the  year  while  a 
high  angle  means  maritime  air  masses  (those  forming  over  large  bodies  of 
water)  dominate.  Perfect  continental ity  and  perfect  oceanity  are  repre¬ 
sented  by  vertical  and  horizontal  lines,  respectively.  His  formula  for 
continental ity  is: 

C  =  LcosA  (7) 

where  L  is  the  comparative  length  of  the  long  axis  and  A  is  the  axis' 
angular  deviation  from  the  vertical. 

By  introducing  precipitation,  Oliver  unnecessarily  complicated  the 
elucidation  of  continental ity.  Annual  precipitation  amount  should  not 
be  used  as  a  measure  of  continental ity,  because  the  occurrence  of  pre¬ 
cipitation  is  influenced  by  more  important  factors  such  as  orography  and 
transitory  circulation  systems. 
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Berg  (1944)  formulated  an  index  based  on  the  relative  frequencies 
of  continental  and  maritime  air  masses: 

K  =  c/(c  +  m)  (8) 

where  c  and  m  are  the  percentage  frequencies  of  continental  and  maritime 
air  masses,  respectively.  Quite  often,  however,  it  is  difficult  to 
distinguish  air  mass  types.  They  often  have  a  combination  of  conti¬ 
nental  and  marine  characteristics,  and  such  a  distinction  is  both  dif¬ 
ficult  and  arbitrary. 

b.  Evaluation  of  existing  indexes 

• 

Contemporary  representations  of  continental ity,  such  as  are  found 
in  textbooks  of  meteorology  and  climatology,  use  indexes  based  on  the 
annual  range  of  temperature  as  noted  above.  The  fundamental  contributor 
to  annual  range,  however,  is  not  continental ity  but  the  variation  of 
insolation  over  the  year.  To  eliminate  this  influence,  these  indexes 
incorporate  a  corrective  term  as  noted  (aS,  sine  0  or  sine  (0  +  10)), 
which,  when  used  as  the  denominator  in  the  expression  for  continental i- 
ty,  presumable  makes  the  resultant  index  essentially  independent  of 
latitude  so  that  only  the  effects  of  continental ity  remain.  As  Conrad 
and  Poliak  (1962)  note,  the  annual  temperature  range  must  be  reduced  to 
equality  at  all  latitudes. 

This  compensation,  however,  is  not  adequate.  Assuming  a  reasonable 
average  atmospheric  transmission  coefficient  of  0.7  and  taking  the  dif¬ 
ference  (AS)  between  summer  (June  +  July)  and  winter  (December  + 

January)  insolation,  then  dividing  by  sin0  or  sin(0  +  10),  yields  values 


that  vary  considerably  with  latitude  (Table  1).  These  values  should  be 
nearly  equal  if  incorporating  latitude  properly  corrects  for  the  lati¬ 
tudinal  differences  in  insolation.  Employing  a  sine  function  implies 
that  AS  is  a  maximum  at  90  degrees,  but  the  highest  value  occurs  at 
approximately  55  degrees.  A  better  correction  could  be  to  divide  annual 
temperature  range  by  the  AS  value  corresponding  to  a  location's  lati¬ 
tude. 

Another  deficiency  of  previous  studies,  in  this  case  of  those  using 
the  lag  of  temperature  behind  radiation  to  indicate  continental ity,  is 
the  implication  (and  perhaps  the  assumption)  that  these  lags  are  equal 
in  length  throughout  the  year.  In  fact,  this  lag  varies.  Annual 
temperature  curves  for  stations  outside  the  tropics  are,  in  general,  not 
symmetrical;  the  lag  of  temperature  after  insolation  is  not  the  same  for 
summer  maximum  and  winter  minimum  (Driscoll,  1984).  Thus,  the  best-fit 
method  of  Fourier  analysis  which  Prescott  and  Collins  (1951)  attempted, 
which  predicated  lag  based  on  the  first  harmonic  only,  is  misleading. 
Continental ity  can  thus  be  seen  as  a  dynamic,  not  static,  climatological 
indicator. 

A  more  appropriate  index  of  continental ity  based  on  lag  should 
therefore  be  possible.  Lags  for  both  winter  and  summer  can  be  construc¬ 
ted  and  mapped.  The  data  needed  are  simply  the  times,  to  the  nearest 
day,  when  mean  daily  temperatures  are  a  maximum  and  minimum  over  the 
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Table  1 

+  10)). 


Latitudinal  variation  of  AS,  AS/sin  0,  and  AS/(sin  (0 


AS  ( Jun  +  Jul ) 
-  (Dec  +  Jan) 


Dividing  AS 
by  Sin  (0) 


Dividing  AS  by 
Sin  (0  +  10) 


90 

791 

791 

803 

85 

793 

796 

796 

80 

802 

814 

802 

75 

823 

852 

826 

70 

865 

921 

878 

65 

915 

1010 

947 

KCv 


3.  OBJECTIVES 


The  objectives  of  this  research  are: 

(1)  To  develop  an  index  of  continental ity  based  on  annual  temper¬ 
ature  range  which  compensates  for  the  annual  variation  of  insolation 
more  satisfactorily  than  existing  indexes; 

(2)  To  develop  indexes  of  continental ity  based  on  the  lag  of 
temperature  behind  insolation  in  summer  and  winter; 

(3)  To  obtain  suitable  temperature  data  for  an  appropriate  number 
of  stations  in  North  America  and  map  the  indexes  developed  in  (1)  and 
(2)  for  this  continent; 

(4)  To  obtain  temperature  data  sufficient  for  a  mesoscale  (region¬ 
al)  analysis  of  continental ity  in  the  Rocky  Mountains  of  the  United 
States  and  southern  Canada,  and  the  area  from  the  Great  Lakes  south¬ 
eastward  to  the  Atlantic  Ocean,  and  therewith  determine  the  influence  of 
elevation  and  the  proximity  of  large  water  bodies  on  continental ity;  and 

(5)  To  compare  and  contrast  the  resultant  maps  of  continental ity, 
attempt  a  physical  (meteorological)  explanation  for  them,  and  evaluate 
each  as  an  appropriate  measure  of  continental ity. 


4.  DATA  AND  PROCEDURE 


For  the  purposes  of  this  study  an  appropriate  distribution  of  sta¬ 
tions  throughout  the  three  political  entities  in  North  America  is  re¬ 
quired.  The  greatest  number  of  stations  with  long-term  temperature 
records  is  in  the  United  States.  Some  of  these  U.S.  stations  were 
deleted  to  ensure  a  more  equitable  distribution  of  stations  throughout 
North  America.  There  were  77  U.S.  stations,  42  Canadian  stations,  and 
9  Mexican  stations  for  a  total  of  128  (Fig.  2)  for  the  macroscale 
(continental)  analysis.  For  the  mesoscale  analysis  of  continental ity  in 
the  Rocky  Mountains,  90  stations  were  selected  (Fig.  3).  About  half  of 
the  stations  were  below  1500  meters  elevation.  One  hundred  nineteen 
stations  were  used  to  analyze  the  Great  Lakes-Appalachia  region  (Fig. 

4). 

There  is  a  question  of  whether  to  use  the  mean  monthly  maximum, 
average,  or  minimum  temperatures  for  each  station.  A  preliminary 
analysis  showed  that  there  is  practically  no  difference  between  using 
maximum,  minimum,  or  average  temperature;  the  aT  values  required  to 
achieve  objective  (1),  and  the  times  of  the  year  at  which  temperature 
maximizes  and  minimizes  (objective  2),  are  virtually  identical  for  all 
three.  Because  mean  monthly  maximum  and  minimum  temperature  records  are 
not  readily  available  for  all  stations  in  North  America,  monthly  temper¬ 
ature  means  were  used  to  approximate  daily  mean  temperatures.  Informa¬ 
tion  for  each  station  consisted  of  30-year  values  of  the  mean  daily 
average  temperature  for  each  month,  as  well  as  latitude  and  elevation. 

Data  should  be  from  coincident  time  periods  for  all  three  countries 
in  North  America.  The  normal  period  1941-1970  was  used  for  Canada 


because  it  was  the  only  period  for  which  temperature  information  was 
available.  The  Mexican  temperature  normals  utilized  were  based  on  the 
period  1931-1960.  Temperature  records  from  the  periods  1941-1970  and 
1951-1980  were  available  for  the  United  States.  The  period  1941-1970 
was  chosen  for  U.S.  stations  to  provide  the  most  overlap  (20  years) 
among  the  three  countries,  thereby  keeping  the  data  as  homogeneous  as 
possible. 

The  first  objective  is  to  develop  an  index  of  continental ity  that 
compensates  for  the  annual  variation  of  radiation  better  than  existing 
indexes.  The  annual  temperature  range  for  each  station  was  divided  by  a 
seasonal  difference  in  radiation  received  (AS)  that  corresponds  with  the 
whole  degree  of  latitude  nearest  that  station.  Values  of  AS  were  calcu¬ 
lated  by  subtracting  the  sum  of  December  and  January  radiation  from  the 
sum  of  June  and  July  radiation.  The  difference  between  January  and  July 
radiation  could  have  been  used,  as  Brunt  (1924)  did,  but  the  mean 
temperature  of  a  particular  month  is  determined,  in  part,  by  the  radia¬ 
tion  received  during  the  preceding,  as  well  as  current,  month.  Since 
most  maximum  and  minimum  temperatures  in  North  America  occur  in  July  and 
January,  respectively,  radiation  received  during  the  months  of  December 
and  January,  and  June  and  July,  was  used  to  calculate  AS.  Regardless  of 

how  as  is  calculated,  the  latitudes  of  maximum  AS  are  almost  identical 

-2  -1 

(Fig.  5).  Since,  when  AS  is  in  units  of  W  m  day  and  range  is  in 
degrees  Celsius,  the  resultant  values  are  small,  they  were  multiplied  by 
a  scaling  factor  of  1000  to  yield  integers.  These  numbers  were  then 
plotted  and  an  isopleth  analysis  performed. 

An  alternative  approach  to  developing  a  better  continental ity  index 
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Fig.  5.  (Jun  +  Jul )  -  (Dec  +  Jan)  radiation  versus  (Jun  +  Jul  + 
Aug)  -  (Dec  +  Jan  +  Feb)  radiation  versus  (Jun  -  Dec)  radiation. 


based  on  annual  range  of  temperature  is  that  of  regressing  the  annual 
temperature  range  on  latitude  and  isolating  the  residuals  from  the  re¬ 
gression  line.  Residuals  were  plotted  and  an  isopleth  analysis  perform¬ 
ed.  Regression  results  are  relative.  Positive  values  can  be  considered 
to  represent  continental  stations  while  negative  values  indicate 
oceanic,  or  less  continental,  climates. 

To  calculate  the  summer  and  winter  lag  of  temperature  behind  radia¬ 
tion,  a  method  of  determining  the  dates  of  maximum  and  minimum  temper¬ 
ature  had  to  be  found.  It  had  to  be  one  that  can  approximate  daily 
temperature  values  from  monthly  means  because  daily  means  are  not  avail¬ 
able  for  North  American  stations  during  the  period  1941-1970. 

The  technique  known  as  cubic  spline  interpolation  was  employed. 
References  to  spline  functions  are  found  in  the  mathematical,  but  not 
the  meteorological,  literature.  A  smooth  curve  is  fitted  through  a 
given  set  of  data  points  on  an  X-Y  plot.  From  that  curve,  values  of  Y 
lying  between  the  given  data  points  can  be  calculated.  In  our  case  the 
X-coordinate  is  time  in  days,  the  Y-coordinate  is  temperature  (°C),  and 
the  interpolating  interval  is  one  day. 

Each  interpolation  interval  is  connected  by  a  third  order,  twice 
differentiable  polynomial.  Second  derivatives  are  then  calculated  and 
joined  together  to  form  an  interpolating  curve  called  a  cubic  spline 
function.  Details  on  cubic  splines  are  found  in  Appendix  A.  Mean  daily 
temperatures  were  generated  for  the  periods  December-January-February 
and  June-July-August,  and  dates  of  maximum  and  minimum  temperature  were 
determined  by  inspection. 

With  the  spline  technique,  a  minimum  of  three  initial  data  points 


is  required  to  determine  dates  of  maximum  or  minimum  temperature.  But 
three  initial  data  points  produced  dates  that  deviated  from  the  "true" 
dates  of  maximum  or  minimum  temperature  by  a  week  or  more.  The  true 
dates  of  maximum  and  minimum  temperature  are  those  dates  on  which  the 
highest  and  lowest  mean  daily  temperatures  usually  occur,  based  on  long¬ 
term  climatological  normal  periods  (e.g.,  30  years).  To  determine  the 
true  dates  precisely  would  entail  using  the  actual  recorded  daily  maxi¬ 
mum  and  minimum  temperatures  and  averaging  them  over  30  years.  Calcula¬ 
ting  true  dates  for  every  station  used  in  this  investigation  would  have 
taken  a  considerable  length  of  time.  In  addition,  daily  temperature 
records  were  not  available  for  all  of  the  stations,  in  particular  those 
in  Canada  and  Mexico.  The  actual  dates  of  maximum  and  minimum  temper¬ 
atures  were  therefore  estimated. 

To  estimate  these  true  dates  using  only  three  data  points,  smooth 
curves  were  drawn  manually  through  the  three  highest  and  three  lowest 
mean  monthly  temperatures  of  a  few  selected  stations,  and  the  dates  of 
maximum  and  minimum  temperatures  determined  by  inspection.  Assuming 
that  the  mean  temperature  of  any  month  occurs  on  the  15th  day,  and  if 
the  preceding  and  the  following  months  have  the  same  mean  temperature, 
then  the  maximum  or  minimum  point  of  the  curve  occurs  on  the  30th  day  of 
the  60-day  interval.  The  spline  function,  however,  places  the  maximum 
(minimum)  temperature  eight  days  later  than  it  should  be.  For  example, 
if  we  assume  that  June  and  August  have  the  same  mean  temperature  (less 
than  July's),  then  the  maximum  temperature  of  the  year  occurs  on  or 
about  July  15,  as  opposed  to  the  spline  date  of  July  23.  After  some 
experimentation,  it  was  found  that  using  five  data  points  (mean  monthly 


temperatures)  instead  of  three  yielded  much  better  results.  Dates  of 
maxima  and  minima  were  consistently  closer  to  true  values,  differing 
non-systematically  by  3-4  days  at  the  most.  It  was  decided,  therefore, 
to  use  the  mean  temperature  of  the  warmest  (coldest)  month  of  the  year 
plus  the  mean  temperatures  of  the  two  preceding  and  the  two  succeeding 
months,  as  input  to  the  cubic  spline  procedure.  If,  for  example,  the 
coldest  month  is  January,  then  the  mean  temperatures  of  November, 
December,  January,  February,  and  March  are  used  to  calculate  the  date  of 
the  minimum  temperature  of  the  year. 

Once  the  dates  of  maximum  and  minimum  temperature  were  known, 
differences  in  days  were  taken  between  those  dates  and  the  dates  of  the 
two  solstices  to  find  the  summer  and  winter  lag  of  temperature  behind 
radiation.  The  differences  were  plotted  and  analyzed.  Lag  patterns 
were  also  compared  with  the  patterns  of  other  continental ity  indexes. 
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5.  RESULTS 


If  we  envision  a  hypothetical  circular,  flat  land  mass  totally 
surrounded  by  water,  it  is  reasonable  to  assume  that  the  highest  conti - 
nentality  is  situated  at  the  exact  center  of  that  land  mass.  Prevailing 
wind  currents  such  as  the  westerlies  would  shift  the  location  of  the 
maximum  continental ity  in  the  direction  of  the  wind  flow  (downstream). 
Topographic  features  affect  the  wind  flow  and  consequently  play  a  role 
in  determining  a  particular  location's  index  of  continental ity. 

Considering  the  North  American  continent  as  a  whole  and  taking  into 
account  the  prevailing  westerly  wind  flow,  the  topography,  and  the  fact 
that  all  of  northern  Canada  is  an  archipelago,  the  theoretical  center  of 
conti nentality  would  likely  be  located  somewhere  along  the  U.S. -Canadian 
border  or  just  to  the  north  of  it,  between  50-55°N  latitude  and  95-100°W 
longitude.  In  the  absence  of  the  Rocky  Mountains,  the  center  would  be 
shifted  eastward  to  an  area  situated  north  of  the  Great  Lakes  and  south 
of  Hudson  Bay. 

An  isopleth  analysis  of  continental ity  index  values  derived  by 
dividing  annual  temperature  range  by  annual  variation  of  insolation  for 
128  stations  in  North  America  is  shown  in  Fig.  6.  The  patterns  exhibit 
some  similarity  to  Conrad's  map  of  continental ity  for  North  America 
(Fig.  1).  An  "X"  on  the  map  means  a  maximum  value  while  an  "N"  signi¬ 
fies  a  minimum  value.  Index  values  are  highest  in  the  northern  interior 
of  the  continent,  with  a  ridge  running  down  the  center  from  north 
central  Canada,  through  the  Great  Plains  of  the  United  States,  into  the 
interior  of  Mexico;  and  lowest  along  the  coastlines,  with  lower  values 
downwind  (to  the  east)  of  the  Great  Lakes  and  Hudson  Bay.  The  tight 
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gradient  along  the  west  coast  reflects  the  blocking  action  of  the  coas¬ 
tal  ranges  against  the  intrusion  of  maritime  air--although  to  a  lesser 
extent  in  the  Puget  Sound  and  Columbia  River  Valley  areas,  where  topo¬ 
graphic  breaks  permit  a  farther  eastward  penetration  of  this  air. 

The  maximum  continental ity  index  values  noted  above  (A/AS)  are  at 
65-70°N,  which  is  also  where  Conrad's  maximum  continental ity  of  North 
America  occurs.  An  index  of  continental ity  that  compensates  for  the 
annual  variation  of  insolation  by  dividing  annual  temperature  range  by 
AS  therefore  does  not  show  the  highest  degree  of  continental ity  where  it 
is  hypothesized  to  be. 

The  next  approach  was  to  regress  annual  temperature  range  on  1 ati - 

• 

tude  for  the  same  128  stations  in  North  America  (Fig.  7).  An  isopleth 
analysis  of  the  residuals  (Fig.  8)  is  similar  to  Conrad's  North  American 
continental ity  index  patterns  (Fig.  1)  and  to  Fig.  6  (A/AS).  In  Fig.  8, 
positive  residuals  indicate  continental  climates  and  negative  residuals 
represent  locations  with  maritime,  or  less  continental,  climates. 
Particularly  evident  in  all  three  of  Figures  6,  7,  and  8  is  the  influ¬ 
ence  of  Hudson  Bay  and  the  Great  Lakes.  Index  values  are  higher  to  the 
west  and  lower  on  the  downwind  (east)  side  of  both  these  water  bodies. 
There  is,  however,  one  notable  exception  to  the  similarity.  While 
Conrad's  maximum  continental ity  and  that  of  A/ AS  are  located  at  65-70°N, 
the  largest  positive  residuals  in  Fig.  8  are  located  at  50-55°N  along 
the  100th  meridan,  which  is  where  the  highest  degree  of  continental ity 
should  theoretically  be. 

Figures  9  and  10  show  the  summer  and  winter  lags  (in  days)  of  the 
maximum  and  minimum  temperature  behind  radiation,  respectively.  If  lag 


Fig.  7.  Annual  temperature  range  (A)  as  a  function  of  latitude  (0)  for  128  stations  In  North 
America,  and  the  resulting  regression  line. 
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were  to  indicate  degree  of  continental ity  correctly,  it  ought  to  leng- 
then  inversely  with  continental ity  (e.g.,  continental  climates  would 
have  the  shortest  lags).  The  only  areas  that  show  the  expected  results 
are  along  the  west  coast  and  in  Florida,  which  have  relatively  long 
summer  lags. 

It  is  quite  clear  that  the  summer  and  winter  lag  patterns  for  North 
America  otherwise  bear  little,  if  any,  resemblance  to  the  corresponding 
conventional  continental ity  index  patterns.  That  is  not  to  say,  how¬ 
ever,  that  conventional  indexes  are  wrong.  But  they  could  be,  in  light 
of  their  deficiencies.  There  are  evidently  factors  other  than  conti¬ 
nental  ity  present.  Some  of  the  summer  and  winter  lag  patterns  can  be 
interpreted  with  the  aid  of  surface  streamline  analyses  of  the  mean  re¬ 
sultant  wind  fields  for  the  summer  and  winter  months  (Wendland  and 
Bryson,  1981;  and  Bryson  and  Hare,  1974).  Topography  and  nearby  bodies 
of  water  must  also  be  considered. 

Summer  lag  is  not  well  differentiated  (Fig.  9).  About  95%  of  the 
continent  has  a  25-35  day  lag.  Exceptions  are  the  west  coast,  the  shore 
along  the  Gulf  of  Mexico,  and  Alaska.  In  the  summer,  the  North  Pacific 
anticyclone  advects  maritime  air  all  along  the  west  coast  of  the  North 
American  continent,  where  the  stations  do  not  reach  maximum  temperature 
until  September  in  some  places.  To  the  east  of  the  coastal  ranges  the 
lag  steadily  decreases.  Clockwise  circulation  around  the  North  Atlantic 
anticyclone  brings  in  moist  air  from  the  Gulf  of  Mexico  into  the  central 
United  States,  causing  the  35-day  isoline  along  the  Gulf  Coast  to 
protrude  northward  toward  the  Great  Plains.  The  moderating  influence  of 
the  Great  Lakes  is  shown  by  the  bulge  or  tongue-like  pattern  over  those 


large  bodies  of  water,  with  the  tongue  pointed  downwind  toward  New 
England.  Some  surface  streamlines  of  tropical  air  run  parallel  to  the 
east  coast  all  the  way  up  to  the  maritime  provinces  (Fig.  11).  As  a  re¬ 
sult,  dates  of  maximum  temperature  along  the  eastern  seaboard  are  some¬ 
what  delayed  relative  to  more  inland  stations,  but  not  delayed  as  much 
as  west  coast  stations. 

A  possible  reason  for  the  short  lag  times  in  the  Yukon  Territory 
and  the  interior  of  Alaska  is  the  presence  of  a  warm  belt  (temperature 
ridge)  at  850  millibars  that  forms  in  April  and  then  vanishes  in 
September  (Bryson  and  Hare,  1974).  The  warm  belt  is  most  likely  the  re¬ 
sult  of  intense  radiational  heating  of  the  plateau  and  mountain  sur¬ 
faces.  As  a  consequence  of  this  heating,  maximum  temperatures  occur 
early  with  respect  to  the  summer  solstice. 

The  winter  lag  pattern  (Fig.  10)  is  more  varied  than  the  summer 
pattern.  Lag  generally  increases  with  latitude,  from  25  days  in  central 
Mexico  to  55  days  in  the  Canadian  archipalego.  Especially  well  marked 
is  the  gradient  in  Alaska.  The  tight  gradient  is  most  likely  due  to  the 
interaction  between  Alaska's  topography,  the  Pacific  Ocean,  and  the 
nearby  polar  ice  pack. 

Several  exceptions  to  the  northward  increase  in  winter  lag  are 
present.  Minimum  values  from  northern  Florida  to  Virginia  could  be 
associated  with  an  anticyclonic  eddy  (region  of  divergence)  in  the 
southeastern  United  States  (Fig.  12).  The  eddy  is  usually  the  remnant 
of  an  outbreak  of  Arctic  or  polar  air,  but  occasionally  the  air  is  of 
Pacific  origin.  Whatever  the  source,  the  air  that  settles  over  that 
section  of  the  country  is  more  continental  (drier  and  colder)  than  the 


surrounding  air  masses,  causing  minimum  temperatures  to  occur  earlier 
than  would  otherwise  be  expected.  The  aridity  of  the  southwestern 
United  States  could  be  a  factor  in  explaining  the  short  lags  in  that 
region.  There  is  less  cloud  cover  to  absorb  the  outgoing  terrestrial 
radiation  and  consequently  the  ground  and  the  air  immediately  above  it 
cool  faster  than  if  the  air  was  more  humid.  Hudson  Bay  begins  to 
freeze  in  mid-December,  but  is  not  completely  frozen  over  until  late 
January.  The  delay  of  the  freeze-over  is  long  enough  to  warm  the  pre¬ 
vailing  northwesterly  Arctic  airstreams  that  pass  overhead,  thereby  in¬ 
creasing  the  lag  to  late  January  for  stations  east  and  southeast  of 
Hudson  Bay.  From  late  January  to  late  June  the  bay  remains  frozen  over, 
making  it  essentially  a  land  surface.  Air  masses  passing  overhead  are 
not  moderated,  as  they  would  be  at  other  times.  Occurrence  of  minimum 
temperatures  are  also  delayed  by  a  few  days  in  the  vicinity  of  the  Great 
Lakes. 

From  a  mesoscale  analysis  of  stations  in  the  Rocky  Mountain  region 
of  the  United  States  and  southern  Canada,  there  is  no  apparent  influence 
of  elevation  on  the  lag  of  temperature  behind  radiation,  for  both  the 
summer  and  winter  cases  (Figs.  13  and  14).  The  other  region  to  be  ex¬ 
amined  more  closely  extends  from  the  Great  Lakes  southeastward  to  the 
Atlantic  Ocean  (Figs.  15  and  16).  Areas  to  the  east  and  southeast  of 
the  Great  Lakes  have  lags  longer  than  would  be  expected  in  the  lakes' 
absence,  with  the  influence  of  those  lakes  being  much  more  pronounced  in 


I 

’  Fig.  15.  Great  Lakes-Appalachia  summer  1 

•  maximum  behind  radiation  maximum. 


6.  SUMMARY 


In  this  research,  an  attempt  has  been  made  to  quantify  conti nenta- 
lity  better  than  has  previously  been  the  case.  First,  conventional  in¬ 
dexes  of  conti nentality  were  examined  and  their  deficiencies  noted. 
Several  approaches  were  then  taken  to  develop  alternative  continental ity 
indexes.  These  approaches  were:  (1)  dividing  annual  temperature  range 
(A)  by  the  summer-winter  difference  in  radiation  receipts  (AS);  (2)  re¬ 
gressing  annual  temperature  range  on  latitude,  and;  (3)  determining  the 
summer  and  winter  lag  of  temperature  behind  radiation.  In  the  first 
approach,  AS  was  used  in  place  of  sin0  and  sin (0  +  10),  the  conventional 
correction  factors,  to  produce  a  continental ity  ratio  which  should  cor¬ 
rect  for  the  differences  with  latitude  of  summer-winter  differences  in 
radiation  receipts.  Linear  regression  was  employed  to  show  the  rela¬ 
tionship  between  latitude  and  annual  temperature  range,  and  to  derive  an 
index  of  continental ity  based  upon  the  deviation  of  annual  temperature 
range  from  the  regression  line.  Summer  and  winter  lag  of  temperature 
behind  radiation  (the  number  of  days  that  the  maximum/minimum  tempera¬ 
ture  occurs  after  summer/winter  solstice)  were  determined  by  the  method 
of  cubic  spline  interpolation. 

The  resultant  maps  of  the  developed  continental ity  indexes  were 
compared  with  each  other  and  contrasted  with  maps  of  conventional  in¬ 
dexes.  Each  of  the  newly  developed  indexes  was  evaluated  as  an  appro¬ 
priate  measure  of  continental ity  and  a  physical  (meteorologica)  explana¬ 
tion  attempted  for  the  patterns. 
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7.  CONCLUSIONS 

Dividing  annual  temperature  range  by  as  produces  a  map  of  conti - 
nentality  that  is  similar  to  maps  of  existing  indexes  of  continental ity. 
A  continental ity  index  with  annual  temperature  range  and  annual  varia¬ 
tion  of  radiation  as  the  parameters  therefore  appears  to  be  no  better 
than  conventional  measures  of  continental ity. 

An  empirical  approach  is  an  index  of  continental ity  derived  by  re¬ 
gressing  annual  temperature  range  on  latitude,  and  then  isolating  the 
residuals  from  the  regression  line.  The  higher  above  the  regression 
line  a  point  is,  the  more  continental  the  location;  the  lower,  the  more 
maritime.  The  region  of  the  highest  residual  values  in  North  America 
comes  very  close  to  the  intuitive  center  of  conti nentality  hypothesized 
previously.  This  index  comes  closest  to  conforming  to  the  hypothesized 
distribution  of  continental ity.  The  regression  technique  takes  into 
account  the  annual  variation  of  radiation  and  apparently  compensates  for 
it  more  adequately  than  existing  indexes  of  continental ity. 

North  American  maps  of  A/ AS  and  regression  residuals  both  have 
similarities  to  maps  of  conventional  indexes  (e.g.,  Conrad).  Isopleth 
analyses  of  the  index  values  reveal  the  influence  of  the  prevailing 
westerly  wind  flow  and  the  topography  as  well  as  large  water  bodies. 

The  A/ AS  patterns  bear  the  closest  resemblance  to  Conrad's  patterns. 
About  the  only  difference  between  the  residual  patterns  and  Conrad's 
patterns  is  that  the  highest  residual  values  are  much  closer  to  the 
center  of  the  continent  than  are  Conrad's  highest  values. 

The  continental  scale  summer  and  winter  lag  patterns  for  North 
America  reveal  to  a  limited  extent  the  influence  of  the  atmospheric  and 


oceanic  general  circulation,  air  masses,  topography,  and  large  bodies  of 
water.  Separate  mesoscale  analysis  of  the  Rocky  Mountain  and  Great 
Lakes-Appalachia  regions  indicate  that  elevation  is  not  a  factor  in  de¬ 
termining  lag.  The  lag  patterns  show  little  resemblance  to  conventional 
index  patterns  or  to  the  patterns  of  A/AS  and  regression  residuals,  and 
thus  other  factors  may  be  present  that  were  not  addressed  in  this  study. 
An  index  based  on  the  lag  of  temperature  behind  radiation  fails  to  place 
the  maximum  continental ity  of  North  America  where  it  should  theore¬ 
tically  be. 

This  research  has  made  it  apparent  that  continental ity  is  relative, 
and  is  a  dynamic,  rather  than  static,  phenomenon  in  that  we  have  shown 
it  varies  with  time  (season)  and  definition.  Of  all  the  indexes  pro¬ 
posed  in  this  study,  the  one  developed  by  regressing  annual  temperature 
range  on  latitude  is  offered  as  an  alternative  to  conventional  indexes 
of  continental ity. 


8.  RECOMMENDATIONS 


There  appears  to  be  an  upper  limit  (envelope)  to  annual  temperature 
range  as  a  function  of  latitude  (Fig.  7).  The  line  containing  the  maxi¬ 
mum  values  of  annual  temperature  range  is  curved,  suggesting  that 
second-order  or  higher  regression  may  be  a  feasible  approach  to  develop¬ 
ing  an  index  of  continental ity.  A  quadratic  solution  possibly  could 
provide  a  better  fit  through  all  the  points  than  can  a  straight  line. 
Examining  other  areas  of  the  world  would  be  helpful  in  establishing  the 
validity  and  usefulness  of  the  continental ity  indexes  developed  in  this 
study,  compared  to  existing  indexes  of  continental ity.  An  investigation 
of  the  Eurasian  land  mass  would  be  particularly  interesting  and  reveal¬ 
ing,  in  light  of  the  enormous  size  of  that  area.  Other  avenues  of  re¬ 
search  are  to  look  at  temperature  variability  and  to  find  more  satis¬ 
factory  explanations  for  the  lag  patterns. 
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CUBIC  SPLINE  INTERPOLATION 


The  most  basic  curve-fitting  problem  is  fitting  a  smooth  function 
through  a  set  of  X,Y  pairs.  There  is  an  infinite  set  of  solutions  to 
this  problem,  two  of  which  are  polynomials  and  cubic  splines  (Cline, 
1975).  Given  n  data  points,  a  polynomial  of  degree  n-1  can  be  passed 
through  them.  One  of  the  difficulties  with  conventional  polynomial 
interpolation  is  that  the  polynomial  may  oscillate  while  the  function 
itself  varies  smoothly.  Another  disadvantage  is  the  extreme  dependency 
of  the  entire  curve  on  each  data  point.  A  slight  movement  of  a  point  at 
one  end  can  drastically  affect  the  curve  all  the  way  to  the  other  end. 

A  way  of  overcoming  these  problems  is  to  use  piecewise  low-order  inter¬ 
polating  polynomials  on  subintervals  of  the  given  interval. 

Such  polynomials  are  called  spline  functions,  introduced  by  I.  J. 
Schoenberg  in  1946.  First  order  spline  functions  are  the  simplest. 

They  are  piecewise  linear  and  not  very  smooth,  but  very  useful  if  the 
spacing  between  nodes  or  knots  (the  given  X  values)  is  small.  Second 
order  splines  lack  symmetry  in  their  determination  with  relation  to 
endpoints  of  the  interval.  Third  order  or  cubic  splines  are  the  most 
widely  used.  An  interpolating  cubic  between  two  adjacent  points  can  be 
denoted  as: 

F,(X)  =  a0  +  a,X  *■  a/  *  a3X3. 

A  smooth  spline  fit  is  achieved  by  connecting  each  pair  of  adjacent 
points  with  a  third  degree  polynomial,  then  matching  up  the  sections  so 
that  the  first  and  second  derivatives  (slope  and  curvature)  are 


continuous  at  each  point  X  and  over  the  entire  region  Xg  <^.  lXn>  The 
approximating  function  over  the  whole  region  is  a  cubic  spline  function 
g(X).  Because  the  second  derivative  of  a  cubic  is  a  straight  line,  the 
second  derivative  of  g (X )  varies  linearly  over  each  interval.  The 
second  derivative  at  any  point  X  is  given  by 


g"(x) 


g"(xi) 


(A.l) 


where  X^  IX  -  Xi+]  (Hornbeck,  1975).  Integrating  twice  and  applying 
the  conditions  that  g(X^)  =  f(X^)  and  g(X^+^)  =  ffX^),  it  is  found 
that  for  the  interval  between  X..  and  X^, , 


9(X)  =  F.(X) 


g"^) 
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AXi(Xi+l 


g,,(xi+l} 


(X  -  X.) 


vir  AXi  (X  -  xi> 


+  f(X.) 


+  f(X, 


i+1 


(A. 2) 


where  AXi  =  Xi+1  -  X^.  Equation  A. 2  determines  the  interpolating  cubics 
over  each  interval  for  i  =  0,  1,  .  .  .  ,  n  -  1.  The  second  derivatives 
are  found  by  using  the  derivative  matching  conditions  F 1  (X^ )  =  ( X^ ) 

and  F V ( X ^ )  =  fV_i (X^ )  to  Eq.  A. 2.  These  conditions  mean  that  g ’ ( Xi )  or 
g"(X.)  are  the  same  when  X  is  approached  from  either  side.  For  i  =  1, 

2,  .  .  .  ,  n  -  1  the  conditions  yield  a  set  of  linear  simutaneous  equa¬ 


tions  of  the  form 
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*  tUs-(x,. 


=  6  - 


f(xi+1)  -  f(xi)  f(X.)  -  f(xi_1) 


(Ax.y 


w^TW~^r 


(A. 3) 


If  the  Xi  are  evenly  spaced,  then  Eq.  A. 3  simplifies  to 

+  [4]g"(x.)  +  [i]g"(x.+1 ) 


=  6 


f(Xi+1)  -  2f (X - )  +  f(X1_1 ) 


(AX.)' 


(A. 4) 


There  are  n  -  1  equations  in  n  +  1  unknowns  g"(XQ),  g" ( X^ ) .  .... 
g"(Xn).  Two  other  required  equations  are  derived  by  specifying  the 
boundary  conditions  at  the  endpoints  g"(X0)  and  g"(Xn).  A  natural  cubic 
spline  is  obtained  if  we  let  9"(Xq)  =  9"(Xn)  =  0.  The  complete  set  of 
equations  is  solved  for  g"(X^),  i  =  1 ,  2,  .  .  .  ,  n  -  1 .  Substituting 
g"(X.j )  back  into  Eq.  A. 2  yields  the  desired  interpolated  values  of  Y. 

Five  data  points  provide  greater  resolution.  Since  the  desired 
value  falls  somewhere  between  the  second  and  fourth  month  (day  30  to  day 
60),  adding  the  two  outlying  months  does  not  significantly  affect  the 
maximum/minimum  dates.  Because  of  the  endpoints,  the  outer  two  inter¬ 
vals  (day  0  to  day  30,  and  day  90  to  day  120)  are  less  reliable  than  the 
two  inner  intervals  (days  30  through  90). 

Program  input  was  five  X,Y  pairs  (monthly  temperatures)  and  one  day 
as  the  interpolation  interval.  Output  consisted  of  120  daily  tempera¬ 


tures  . 


There  is  an  error  term  associated  with  the  spline-derived  tempera¬ 
ture  values  which  could  affect  the  lag  results.  The  spline  value  of  the 
mean  temperature  of  any  month  (sum  of  the  daily  values  divided  by  the 
number  of  days  in  that  month)  should  equal  the  actual  mean  temperature. 
But  as  shown  in  Table  A1 ,  the  spline  and  actual  values  are  not  identi¬ 
cal,  as  calculated  from  a  sample  of  ten  stations.  It  is  interesting  to 
note  that  the  differences  between  the  actual  and  spline  values  are  all 
positive  for  the  summer  months  while  the  differences  are  all  negative 
for  the  winter  months.  In  addition,  the  more  continental  a  location  is 
(the  greater  the  amplitude  of  the  annual  temperature  curve),  the  larger 
is  the  difference. 

If  the  actual  maximum  and  minimum  temperatures  occur  on  the  same 
dates  as  the  spline-derived  maximum  and  minimum  temperatures,  then  the 
error  term  would  not  change  the  lag  values  at  all.  On  the  other  hand, 
if  the  dates  of  occurrence  are  not  identical,  the  effect  on  the  lag 
results  could  be  significant.  Since  the  differences  are  proportional  to 
degree  of  continental ity,  the  lag  values  may  change,  but  not  much.  The 
overall  lag  patterns  would  probably  remain  essentially  the  same. 
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983 

Moosonee,  Ont. 

51 .3 

35.7 

985 

Mould  Bay,  N.W.T. 

76.2 

39.3 

817 

Nitchequon,  Que. 

53.2 

36.5 

990 

Norman  Wells,  N.W.T. 

65.2 

44.8 

915 

Nottingham  Is.,  N.W.T. 

63.1 

30.4 

935 

Port  Hardy,  B.C. 

50.7 

11.4 

985 

Prince  George,  B.C. 

53.9 

26.7 

990 

5' 


Station 

Latitude  (°) 

A  (°C) 

AS  (W  m_l 

McGrath,  AK 

63.0 

37.5 

934 

St.  Paul  Is.,  AK 

57.2 

13.6 

987 

Yakutat,  AK 

59.5 

16.2 

970 

Tucson,  AZ 

32.1 

19.7 

746 

Winslow,  AZ 

35.0 

25.4 

803 

Bishop,  CA 

37.4 

21.9 

835 

Eureka,  CA 

40.8 

5.4 

898 

San  Diego,  CA 

32.7 

9.0 

765 

San  Francisco,  CA 

37.6 

8.7 

854 

Santa  Maria,  CA 

34.9 

6.7 

803 

Denver,  CO 

39.8 

24.0 

884 

Grand  Junction,  CO 

39.1 

28.9 

869 

Dulles  IAP,  D.C. 

39.0 

24.0 

869 

Jacksonville,  FL 

30.4 

14.6 

716 

Mi  ami ,  FL 

25.8 

8.7 

619 

Tallahassee,  FL 

30.4 

15.9 

706 

Tampa,  FL 

28.0 

12.1 

663 

Atlanta,  GA 

33.6 

19.8 

784 

Boise,  ID 

43.6 

25.3 

928 

Moline,  IL 

41.4 

29.4 

904 

Ft.  Wayne,  IN 

41.0 

26.5 

898 

Dodge  City,  KS 

37.8 

26.9 

854 

Topeka,  KS 

39.1 

27.9 

869 

Louisville,  KY 

38.2 

24.2 

854 

New  Orleans,  LA 

30.0 

16.1 

706 

Shreveport,  LA 

32.5 

20.0 

756 

Caribou,  ME 

46.9 

30.1 

962 

Portland,  ME 

43.6 

25.8 

934 

Alpena,  MI 

45.1 

26.5 

945 

Int'l  Falls,  MN 

48.6 

35.5 

972 

Saint  Cloud,  MN 

45.6 

34.0 

950 

Jackson,  MS 

32.3 

19.2 

746 

55 
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Station  Latitude  (°) 


A  (°C) 


AS  (W  m-2day-l) 


Station 

Latitude  (°) 

A  (°C) 

AS  (W  nr^dayl ) 

923 


Casper,  WY 


42.9 


26.6 


Conventional  and  proposed  continental ity  index  values  for  the  continental  scale  stations 
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Juan  Manuel  Yee  Fong  was  born  in  Havana,  Cuba  on  1  January  1955. 

He  immigrated  into  the  United  States  in  1962  and  settled  in  Albuquerque, 
NM,  where  he  graduated  from  Highland  High  School  in  1973  and  became  a 
naturalized  U.S.  citizen  in  August  1976. 

In  December  1977,  Juan  received  his  Bachelor  of  Science  degree  in 
Biology  +  Chemistry.  He  was  commissioned  into  the  U.S.  Air  Force  in 
May  1979,  then  was  assigned  to  the  University  of  Texas  at  Austin  from 
June  1979  to  May  1980  for  the  AFIT  (Air  Force  Institute  of  Technology) 
Basic  Meteorology  Program.  His  next  assignment  was  at  Luke  AFB,  AZ  from 
June  1980  to  July  1983,  where  he  served  as  a  forecaster  and  Wing  Weather 
Officer. 

Juan  came  to  Texas  A&M  University  in  August  1983  to  pursue  the 
degree  of  Master  of  Science  in  Meteorology  under  the  sponsorship  of 
AFIT.  Captain  Yee  Fong  is  currently  assigned  to  the  Air  Force  Global 
Weather  Central,  Offutt  AFB,  NE.  His  permanent  mailing  address  is: 

516  Georgia  SE,  Albuquerque,  NM  87108. 
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